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Abstract 

We show that the basic equation of theory of open systems, Gorini- 
Kossakowski-Sudarshan-Lindblad equation, as well as its linear and 
nonlinear generalizations have a natural classical probabilistic inter- 
pretation ~ in the framework of prequantum classical statistical field 
theory. The latter gives an example of a classical probabilistic model 
(with random fields as subquantum variables) reproducing the basic 
probabilistic predictions of quantum mechanics. 

1 Introduction 

Recently interest to probing the limits of quantum mechanics essentially in- 
creased - both quantum contra classical and quantum contra "superquan- 
tum" models (i.e., models deviating from classical physics even more than 
it does "conventional quantum mechanics"), see, e.g., [I]-[8]. On one hand, 
this situation is a consequence of the internal development of physical science. 
On the other hand, it became clear that the projects of quantum computing 
and quantum cryptography have some difficulties in realization and new deep 
foundational analysis is needed. In particular, the role of environment can- 
not be ignored and, instead of unitary dynamics, open systems dynamics (in- 
cluding non-Markovian) may play an important role in quantum computing 
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and other applications. In this paper we present the quantum contra clas- 
sical probing of open system dynamics; in particular, Gorini-Kossakowski- 
Sudarshan-Lindblad (GKSL)equation. The result of our consideration does 
not match with the rather common expectation. We show that, for any semi- 
group of positive (in particular, completely positive) maps, the density oper- 
ator dynamics can be interpreted in the purely classical way. We derive the 
corresponding differential equation for the classical probabilistic dynamics. 
It is very complicated comparing with the original linear operator dynam- 
ics: in general nonlinear and integro-differential. Therefore operationally it 
is natural to work the standard quantum dynamical equation. However, the 
possibility to interpret this equation in the classical probabilistic terms has 
nontrivial foundational consequences. 

The study in this paper has the natural coupling to prequantum clas- 
sical statistical field theory (PCSFT), see, e.g., [9], [3], [7], [8], [10]. The 
latter gives an example of a classical probabilistic model, with random fields 
as subquantum variables, reproducing the basic probabilistic predictions of 
quantum mechanics. The basic idea beyond PCSFT is very simple: to couple 
quantum states, density operators, with the covariance operators of classical 
("prequantum") random fields in complex Hilbert spaces. Of course, the co- 
variance operator does not determine a random field uniquely. However, by 
restricting prequantum fields to Gaussian fields we determine the prequan- 
tum field by its covariance operator. (We consider fields with zero average.) 
Of course, in general the covariance operator of a random field has non-unit 
trace. Therefore on the subquantum level there is no reason to consider the 
trace-preserving dynamics. And we proceed in this paper with semigroups 
of positive maps which in general are not trace preserving. In section H] we 
derive the dynamical equation for density operators corresponding to such 
dynamics of covariance operators. This equation is nonlinear and it has 
quadratic nonlinearity. 

2 Prequantum Classical Statistical Field The- 
ory 

Here we briefly present essentials of PCSFT, see [9], [3], [7], [8], [lO] for 
detailed presentation. To simplify considerations, we will study quantum 
systems with finite dimensional Hilbert spaces. 

Take complex Hilbert space H as space of classical "prequantum" states. 
Consider a probability distribution /j, on H having zero average (it means 
that Jjj{4>, y)dfi{(f)) = for any y E H) and the covariance operator B which 
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is defined by the Hermitian positively defined bilinear form: 

{Byi,y2)= / {yi,(j)){(l),y2)dfx{(j)),yi,y2e H, (1) 
Jh 

By scaling we obtain the operator 

P = Pb^ B/TrB, (2) 

with Trp = 1. Mathematically it has all properties of the density operator. 

By PCSFT a quantum state, a density operator, is simply the symbolic 
representation of the covariance operator of the corresponding prequantum 
(classical) probability distribution. 

In general a probability distribution is not determined by its covariance 
operator. Thus the correspondence PCSFT QM is not one-to-one. How- 
ever, if the class of prequantum probability distributions is restricted to Gaus- 
sian, then this correspondence becomes one-to-one. 

In PCSFT classical variables are defined as functions from (Hilbert) state 
space H to real numbers, / = By PCSFT a quantum observable, a 

Hermitian operator, is simply a symbolic representation of / by means of 
its second derivative (Hessian), f ^ A — i/"(0). This correspondence is 
neither one-to-one. However, by restricting the class of classical variables 
to quadratic forms on Hilbert state space H, fA{<^) — {^(f),(f>), we make 
correspondence PCSFT — )■ QM onc-to-onc. And finally, we present the basic 
equality coupling prequantum (measure-theoretic) and quantum (operator) 
averages: 

if a), = EfA^ [ fA{<f>)dfx{<f>) = TVBA = (TVB) (A)„ (3) 

where 

(A), = TrpA. (4) 

Thus the quantum average {A)p can be obtained as just scaling of the classical 
average. We remark that scaling parameter TrD is, in fact, the dispersion of 
the probability distribution /i : 

a'^EMl'^ [ M\'d^i{<f>)^TrB. (5) 
Jh 

It describes the strength of deviations of a random vector (j) from its average. 
The latter is zero in our model. U H = L2(R'"), the space of square integrable 
functions, : R"^ C, then we can consider as a random field: = 
(p{x,uj), where x e R"* and u is the chance parameter. The average is the 
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field which is equal to zero everywhere. Suppose that is very small, so 
deviations from zero are small. Since 

{A), = ^{U)„ (6) 

the quantum average can be considered as the average of amplified "prequan- 
tum signal". 

Remark. The basic relation of PCSFT coupling quantum states (density 
operators) with covariance operators of prequantum processes, see ([2]), can 
be used the other way around. M. Ohya and N. Watanabe [TT] used it in 
classical signal theory to define information (in the case of Gaussian noisy 
channels with the infinite numeber of degrees of freedom) with the aid of von 
Neumann entropy. 



3 Dynamics of the probability measure cor- 
responding to semigroup of positive maps 

Consider a semigroup of positive maps Ct = e*"^ which is induced by the 
generator given by the superoperator C. Although the general form of the 
generator C is unknown, there are known numerous examples of such gen- 
erators. For examples, the general form of the generator of the completely 
positive semigroup was described in the relation with the GKSL-equation. 
The latter equation gives the most important (for applications) example of 
the dynamics under consideration. We set Bt = CiBq, where Bq is a positive 
operator Hermitian operator. Then Bf is a family of positive operators giving 
the solution of the Cauchy problem: 

dBt 



dt 



CBt, Bt=o = Bq. (7) 



Let us denote by Pt the Gaussian measure with zero mean value and the 
covariance operator Bf. Our aim is to derive the differential equation for 
the map t — )■ Pf. To simplify the derivation, we consider the case of finite 
dimensional Hilbert space. Here Pt has the density 

pt(z, z*) = , ^ e-^B-''''\ (8) 

"^'^ ' v/(2vr)MetBt ^ ' 

By selecting some orthonormal basis (cj), i.e., z = J^j^j^j^^j — (^)^j)) we 
represent the density as 



This density can be represented by means of its Laplace- Fourier transform 
as 

or in the coordinate form: 



J e-(Bti,0+i{i,z)+i{z,Od^dC (10) 



Pt[z,z 



We have 

dpt. 



dt I 

'^ZjOZ.- J 



5^(i?„r(|e.)(e,|))^p,(z,z*), 



dz,;dz. 



where in the last expression we used the scalar product in the space of Hilbert- 
Schmidt operators. We can expand Bt with respect to the orthonormal basis 
(|efc)(em|), 

Bt = y^^{Btek,em)\ek){em\. 

km 

Finally, we have 



dt ^ ' ' 

km;ij 



UkU*^Pt{u,u*)dudu*-^—^Pt{z,z*), (12) 



where 

^km,ij = mek){em\, \ei){ej\) (13) 
By using notation from statistical mechanics we write our equation as 

^(z, z*) = J2 ^km,ij'^-g^g-;Ptiz, z*). (14) 

km;ij * 

We now show that, for any Gaussian solution of the latter equation, its 
covariance operator satisfies the equation ([7]). We have 



z, z*)dzdz* 
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-J (£5(^,0 ( j ZiZ*e'^^^'^+'^'^^^dzdz*y~'^^'^^^^didC 

j / ^'^'^'^""^'^'^dzdz^y^^^'^'^didc 



Thus any dynamics of the form ([7]) in the space of positive operators can 
be imagined as corresponding to the dynamics fll4l) in the space of Gaussian 
measures. The latter is nonhnear and it is integro-differential. Of course, the 
former is essentially simpler than the latter. Therefore one can restrict con- 
siderations to the dynamics of operators and proceed in the phenomenolog- 
ical framework, i.e., without coupling these operators to Gaussian measures 
("distributions of hidden parameters"). 

Example 1. Consider our equation for the probability density in the 
real case and for dimension one. It has the form: 

^(x) = ^x2pi(x). (15) 

(the factor 1/2 is related to consideration of the real case). Consider the 
solution of the form 

3^ 

e 

Then 

dpt, , ^ / 1 \ , , 

lF'^> = t(BF^l^>*<"' 

and ^ 

// / N ^ X 1 



The equation fllSp implies that ^ = ABt, the covariance (in this case simply 
dispersion) evolves linearly. 

Example 2. Consider the dynamics: 

'^ = ABt + BtA\Bt=, = Bo, (16) 

where A is some operator (in general non-Hermitian). Then 

Bt = e^'Boe^*'. (17) 
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Consider now the dynamics of the "state vector" with a random initial con- 
dition 

^{u) = AUu), 0i=o(c^) = M^), (18) 

where 0o has the Gaussian distribution with zero mean value and the covari- 
ance operator Bq. Its solution (pti^j-^) = e^^(po{^) has the covariance operator 
given by f lT7|) . At the same time we know from theory of classical stochastic 
processes that the probability distribution of (pt satisfies the trivial forward 
Kolmogorov equation with zero diffusion. This equation is linear. Thus 
in some special cases nonlinearity and integro-differential structure of the 
equation for the Gaussian probability distribution corresponding to the lin- 
ear dynamics of the covariance operator are redundant. In particular, set 
A = —iH, where H is the hermitian operator. Then the equation (fT6l) is 
nothing else than the von Neumann equation, see also [3]: 

^ = -t[H,Bt],Bt=o = Bo, (19) 

4 Nonlinear dynamics of the density operator 

Our aim was representation of quantum mechanics with the aid of classical 
Gaussian distributions. In general the trace of a covariance operator is not 
equal to one. Moreover, the dynamics considered in the previous section need 
not preserve the trace, i.e., even by starting with a density operator po = Bq, 
we need not obtain a trajectory in the space of quantum states. 

To obtain the density operator corresponding to the covariance operator 
of the prequantum random field, we normalize the covariance operator of a 
Gaussian measure by its trace: 

pt = —4-. (20) 

Then the differential equation ([7]) is transformed to a differential equation 
for p. 



dt TrBt Tr^Bt 
CBt BtTiCBt 

Thus the evolution equation for the density operator corresponding to the 
covariance operator is nonlinear with quadratic nonlinearity: 

^ = Cpt- pt^rCpt. (21) 
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Suppose that the original covariance-dynamics is trace preserving. Than we 
have: 

dTipt dpt 

Thus the last term in (pT]) disappears and dynamics for the covariance and 
density operators coincide and they are both linear: 

^ = Cpt, pt=o = po- (22) 

However, from the prequantum viewpoint there is no reason to consider only 
trace preserving dynamics of covariance operators. Hence, the basis dynam- 
ical equation is the nonlinear equation (12T|) . 



5 Brownian motion and Ornstein-Uhlenbeck 
processes 

Consider the Brownian motion in H. It is determined by its covariance 
operator E (Hermitian and positively defined). The process WY;it,u) = 
(^o(w) + V^w{t, u), where w{t) is the standard Brownian motion, i.e., w{t) ~ 
J\f{0,tl), where / is the unit operator and is a random variable. Suppose 
that the latter is independent from w{t)] denote its covariance operator by 
Bq. Then the covariance operator Bt of wj:{t) satisfies to the equation: 

= Bt=o = Bq. (23) 

The corresponding "density operator" has the form 

_ Bo + J:t 
^* ~ Tr(Bo + St) ■ 

This function satisfies to the ugly differential equation with time dependent 
coefficients. Thus the classical Wiener process is not among processes corre- 
sponding to semi-groups of positive linear maps. 

Now consider the solution of stochastic differential equation 

d(f)t{i^) = At(j)t{u)dt + dwY.{t, w), 0o(w) = ^o('^), 

where the initial condition is again independent from the ws(t). (In particu- 
lar, this type of equations describes the Ornstein-Uhlenbeck process). Then 
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it is easy to show (by using the Ito formula) that its covariance operator Bf 
satisfies the equation: 



- AtBt + BtAl + E, Bt=o = Bq. 



(24) 



dt 



Thus the classical Ornstein-Ulenbeck process is neither in the club of "pre- 
quantum processes." 

By selecting the time independent coefficient At — —1%, where is a 
Hermitian operator, we obtain the dynamics: 



This equation is just inhomogeneous von Neumann equation. (It might be 
that such equations were used in QM, but I do not know such applications.) 

This paper was prepared during the visit to the Center for Quantum 
Bio-Informatics (Tokyo University of Science), February- March 2013. The 
author would fike to thank Massanori Ohya and Noboru Watanabe for fruitful 
discussions and hospitality. 
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